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Abstract 



A recently developed helicity basis for nucleon-nucleon (NN) scattering is 
applied to the deuteron bound state. Here the total spin of the deuteron 
is treated in such a helicity representation. For the bound state, two sets of 
two coupled eigenvalue equations are developed, where the amplitudes depend 
on two and one variable, respectively. Numerical illustrations based on the 
realistic Bonn-B NN potential are given. In addition, an 'operator form' of the 
deuteron wave function is presented, and several momentum dependent spin 
densities are derived and shown, in which the angular dependence is given 
analytically. 
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I. INTRODUCTION 



In a recent manuscript |3J] we developed a three-dimensional approach in momentum 
space for nucleon-nucleon (NN) scattering. The motivation is that for higher energies too 
many partial waves contribute and a direct solution seems more natural and economic. As 
relevant variables momentum vectors appear, specifically their magnitudes and the angles 
between them. The formulation in [|l| is based on a helicity representation with respect to 
the total spin of the two nucleon system. This representation is different from the often 
used helicity basis referring to the individual nucleons |H§. A further important advantage 
of a three dimensional approach is that a sometimes tedious partial wave expansion of a 
complex NN force is no longer needed. Instead one introduces a helicity representation of 
the NN force, which is perfectly adapted to the set of six operators completely describing 
the most general NN force compatible with general invariance principles. Thus, for any NN 
force given in operator form this scheme is applicable. 

The helicity representation developed for NN scattering can also be applied to the bound 
NN system. It may appear unnecessary to extend this particular formulation to study the 
nonrelativistic deuteron, which only contains S and D waves. However, the standard practice 
requires a partial wave representation of the NN force, which we avoid. It is straightforward 
to extend the helicity formulation developed for scattering to investigate the deuteron, cal- 
culate its binding energy, and wave function properties. This is the purpose of the present 
investigation. In addition we study the various wave function properties in momentum space 
in a three dimensional fashion. Our numerical example is based on the Bonn-B potential 
@]. A graphical study of similar character has been carried out in configuration space f| 
based on the AV18 nuclear force ||. In addition we derive and display probability densities 
of various spin configurations for an overall polarized deuteron. These densities are based 
on analytical expressions, which we think are new. 

In Section II we introduce the expansion of the deuteron state into helicity basis states 
as defined in [[L]]. Then we project the deuteron eigenvalue equation on these states. Since 
the deuteron has spin 1 there are three possible values for the helicity projections, namely 
A = 1, 0, —1. Symmetry properties allow to consider only A = 1, 0. Thus we obtain a set of 
two coupled equations in two variables, the magnitude |q| of the relative momentum vector 
and the angle between q and the arbitrarily chosen z-axis. This two dimensional form of the 
deuteron wave function is then connected to the standard partial wave representation. We 
demonstrate that this set of two coupled two dimensional equations can be readily solved 
and display various wave function properties. 

In Section III we derive the deuteron wave function in 'operator form'. In a configuration 
space representation such a form has been given before M. For our purpose an 'operator 
form' is an ideal starting point, since the spin degrees of freedom appear explicitly as spin 
operators, and thus fit perfectly into our helicity formulation. The projection of the wave 
function on the helicity basis leads to deuteron wave function components with an analytical 
angular behavior, which is different from the familiar one. The 'radial' part of the wave 
function satisfies a set of two coupled (A = 1, 0), one dimensional eigenvalue equations in 
|q|. Based on this more analytical insight the connection to the standard 5* and D waves 
forms started in Section II can be finalized. 

There are various possibilities that the two nucleons in the deuteron have a specific 
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orientation of their spins for an overall polarized deuteron. For instance, both nucleons can 
have their spins up, or one nucleon can have its spin up and the other down. In Section 
IV we derive analytic expressions of the corresponding probabilities and display the results. 
This may have applications for electron scattering on the deuteron. Finally we summarize 
in Section V. 

II. FORMULATION I 

A. Deuteron Wave Function in the Helicity Basis 

Let 1^ d ) represent the deuteron state. Here M d is the projection of the total angular 
momentum along a chosen axis, e.g. the z-axis. The deuteron state will now be represented 
in the helicity basis |q; qSA; t) 7 ™ defined in Here q stands for the relative momentum of 
the two nucleons, S for the total spin, A for its projection along q, and t for the total isospin. 
The index ir denotes the parity of the state, and a indicates the state being antisymmetric. 
This results in 

*f d ) = 7 E / rf q |q; <? 1A ; o> la la ( q; $ia ; o| . (2.1) 

Here we took into account that for the deuteron S = 1, t = and the parity is even. The 
general form of the helicity eigenstate is given by JU 

\q;qSA;ty a = (|q> + ^ |-q)) |gSA> \t) , (2.2) 

where r/ n denotes the parity eigenvalue. Using \—qS A) = (— ) s \qS — A) one verifies the 
following properties: 

|q; qSA; t) m = V A- f (|-q> + V, |q» \~qS ~ A) \t) 

= vA-) s \-*-qS-A;ty a . (2.3) 

With the above relations one obtains for the integral in Eq. ( |2.1| ) 

J dq |q; ql - 1; 0) la la (q; ql - 1; 0|*f*) = J dq \ - q; -qll; 0) la la (-q; -gll; 0|*f d ) 

= Jdq |q; gll; 0) la la (q; gll; 0|M>^>. (2.4) 

Hence, Eq. ( p.l|) simplifies to 

K d ) = fdq{^\q;qll;0) la la (q; gll; 0| + i |q; glO; 0) la la (q; glO; 0| 

= / rfq |q; gll; 0) 1 >f rf (q) + \ |q; glO; 0) la y^(q)} , (2.5) 
where we defined 

^f d (q)= la (q;glA;0|vl/f ). (2.6) 
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The azimuthal dependency of the amplitude defined in Eq. ( |2.6| ) can be found as follows. 
The state \q;qSA) is obtained by rotating the state \qz;zSA) from the z-axis into the 
direction of q as 



|q; qSA) = R(q) \qz;zSA) 



(2.7) 



where R(q) = exp —iJ z (fiexp —iJ y 9, and J = L + S is the operator of total angular momen- 
tum. It follows that 



la (q; qSA; t\ *f d ) = la (qz; zSA; t\ e iJ ^ B e iJ ^ 
= e iMd ^ la (qz;zSA;t\e iJ « 



a 



(2.8) 



Thus, we can redefine (p A d (q) such that the azimuthal dependency is factored out 

^f%) = ^(g,0)e^. (2.9) 
This leads to the final expression of the deuteron state in the helicity basis 



a 



rfq { ~ |q; gll; 0) la ^ d (q, 9) + ~ A |q; glO; 0) ia tp%"(q 



la . M d 



JM d <t> 



(2.10) 



The normalization of the wave function components (p^ d (q,9) can be calculated as 



dq'| dq!± la { q ';q'll;0\cp^( q ',9 f ) + ^ la (q'j g'10; 0| <ftf"tf, 9')} 
{i |q; gll; 0) la cpf^q, 9) + \ |q; glO; 0) la #(g, 0)} 
J dc(j d q !^ la (q';g'll;0| q; gll; 0) la ^*(q', 9')^(q,9) 



1 



+-- la <q'; g'10; 0| q; glO; 0) la <^*(g', tf)ttf *{q,9) 
lb 



1 



+- la <q';g'll;0| q; glO; O^^V^'Vo^, 
+ l la (q'; g'10; 0| q; gll; 0) la ^*(g', #'Vi /d (<7, *) 



o 



i 



2vr / dgg 2 / dcosfl ^ <^f d (g,#) 



+ 4 



(2.11) 



Here we used that 

(q'; gW; f | q; gSA; t) m = (1 - i&r(-) 5+ ' 



(«y(q'-q)<JA'A + ^(-) S <y(q' + q)5 



A', -A 



(2.12) 



From the normalization calculated using Eq. (|2.1|) we define a deuteron momentum density 



/0 M< ^qj as 



p Md (q) 



1 



h>\ (q, 



i 



(q) 



4 ^ ^ 4 

1 Jfc,_ m 2 , 1 

4 



1 



<Pi (-q) 



4 

2 1 

+ 4 



(2.13) 
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B. The Deuteron Eigenvalue Equation in the Helicity Basis 



The deuteron state \^ d Id ) satisfies the eigenvalue equation 



{H -E d + V) 



M, 



0. 



(2.14) 



with E d being the deuteron binding energy and V the NN potential. This eigenvalue equation 



is projected on the basis states |q;glA;0) la introduced in the previous Section. Using 
Eq. fl2.10|) for representing 1^ d ) one obtains 

la (q;qlA;0\{H o -E d + V) 



a 



la (q; qlA;0\(H -E d ) 
+ la (q; glA; 0| V J dq' |q'; g'll; 0) la tf'tf, 9') + \ |q'; g'10; 0) lQ #(g', 9')] e* 



m 



-E d )^(q,9)e m ^ 



-~ / dq! la (q; glA; 0| V |q'; g'll; 0) la ifptf, 9>)e 



la ^(q',9')e lM ^ 

la ,J&&( J Q'\JM d (j>' 



-~ I dq' la (q; glA; 0| V |q'; g'10; 0) 1 " ^(g', 9')e 



£ 



(2.15) 



where V^f/(q, q') = na (q; qSA;t\ V |q'; q'SA'; t) na . Thus, the eigenvalue equation for 
deuteron binding energy consists in the helicity basis of a set of two coupled integral equa- 
tions. 



-E d \^(q,9) 



m 



+ J dq'e-^C^') {\v™(q, q>f d (g', 0') + \v£°(q, q')^(q', 9')) = 0, (2.16) 

where the index A corresponds to 1 or 0. Since the wave function components ip^ d (q, 9) have 
no azimuthal dependence, the (//-integral in Eq. (|2.16|) can be carried out independently. 
Defining 



vS Md (q,q',e,e f ) = r#' e -WV A f(q,q') 



(2.17) 



the coupled integral equations are actually only two-dimensional and their final form reads 



E d \^(q,9) 



m 



2 Jo 



dq'q' 2 [ dco S 9'v™> M *(q,q',9,9 f )^(q',9') 



1 pOO rl 

+ ll dq ' ^ J i dC ° s6 ' v ™ & 9 >° d tt' ff ) = °- ( 2 - 18 ) 
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The eigenvalue equation, Eq. ( |2.18|) , is consistent with our treatment of the NN contin- 
uum of Ref. 0, where we derived coupled integral equations in two variables for the NN 
scattering equation. We would like to add the remark that it is not necessary to use any 
specific information about the the deuteron, namely spin, isospin and parity. The set of 
equations, Eq. fl2.16| ), is valid for arbitrary S, 77^, and t. Any calculation based on a realistic 
NN potential will then reveal that a solution of the eigenvalue problem exists only for the 
well known quantum numbers. Therefore, the scheme laid out above automatically provides 
full insight, no a-priori knowledge needs to be employed. At this level the question of total 
angular momentum of that bound state is undetermined and will be considered in Section 
III. 



C. Partial Wave Projection of the Deuteron Wave Function 



In this subsection we want to relate 



^ d d j as given in Eq. (|2.10| ) to the standard partial 
wave representation of the deuteron. The standard representation of the total angular 
momentum basis \q(lS)jm;t) with the normalization 



(q'(l'S')j'm' \q(lS)jm) = Q Hvi&S'S&j'j&m'n 

q'q 



(2.19) 



is used for projecting the wave function ^ d d ). Again, we use the fixed spin and isospin of 



the deuteron. The quantum numbers I, j, m remain arbitrary, and we obtain 



^2.20) 



= \ J dq! (q(ll)jm; 0| q'; g'll; 0) 1 Vf d (q',e')e m ^ 
+\ J dc{ (q(ll)jm; 0| q'; g'10; O) 1 "^^, 9')e iM ^ . 

Recalling the explicit representation of the helicity state ffl, 

\^qSA) = |q> \qSA) 

= £ \q(lS)jm) £ C(lSj; //, m - /i, m)Y^q)e^ m -^d s m ^ A {6), 

Ijm M 

and using Eq. Q2.2|) , the scalar product of partial wave and helicity basis can be worked out 
as 

(q'(lS')jm; t'\ q; qSA; t)* a = {{q'{lS')jm\ q) \qSA) + Vw (q'(lS')jm\ -q) \qSA)) (t>\ t) 

8{q'-q) 



(2.21) 



q'q 



£ C(lSj; ii,m-n, m) (Y^q) + r^(-g)) e^ m ~ ^ d s m _ ^ K {6) 



q'q 



£ C(lSj; ti,m-fJL, m)Y l l{q)e^ m -^d s m _^0) 
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^l + ^(-) jy— ——Spsfae 

M^e*«C{lSr, OAA)Cv(0). (2.22) 



Here we used the relation 



^(?)^ = \/^<oW (2-23) 
together with an addition theorem for Wigner's D-functions 

£ C(ZSj; ^ m - At, ™)<o(0)4- ft A(0) = C(ZSj; OAAR A (0). (2.24) 

Hence, the projection of the deuteron state on the partial wave basis as defined in Eq. (|2.20|) 
is given by 



Ml) = (1 + HO f 1 ^ J** d<p>e-^- M ^' 

1 dcostf' {^C(Zlj; 011)^ ml (^)^ /d (g, f) + \c(llj; OOO)4o(0>^(g, ^ 



(l + (-)') ^(2/ + l)5 mMd 

dcostf' (^(/lj-Oll)^^)^^^) + ^(Zlj-OOOKo^O^M') f (2.25) 



This projection exists only for m = and even /, as enforced by the even deuteron parity, 
and one obtains 



^(q) =2^(21 + 1) 

1 rfcos^' {\c(llj; Oll)4 dl (0>f% ff) + ^C(llj; OOO)4 dO (0')<^ d (g, 6') \ (2.26) 



i 



At this point the remaining properties for the projections ipi(q) can be found by explicit 
calculation. In other words the fact that I has to be and 2, and consequently has to 
be inferred numerically from the solution of Eq. (|2.26 ). Of course an analytical investigation 
can be added once we adopt the deuteron wave function with the well established analytical 
angular behavior (see Section [111 C| ) . 



D. Explicit Solution of the Deuteron Eigenvalue Equation 

The numerical solution of the set of two coupled eigenvalue equations Eq. (|2.18| ) poses 
no specific difficulty We employ an iterative, Lanczo's type technique [HH, which provides 
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both, eigenvalue and eigenvector. Following Ref. the method is modified to avoid un- 
physical solutions corresponding to bound states in the repulsive core region. All numerical 
calculations are carried out with the Bonn-B potential [|J. 

For the discretization of Eq. fl2.18| ) we employ Gaussian grids. The 0' integration over 
the potential, Eq. ( |2.17| ), needs only 10 quadrature points, whereas the cos 9' integration 
requires at least 32 grid points and the q' integration 72 grid points, depending on the 
desired accuracy. The q' integration interval can safely be cut off at 30 fin -1 . Using these 
numerical parameters, we obtain the deuteron binding energy 2.224 MeV. 

In Fig. |TJ the two deuteron wave function components (p^ d (q, 9), A = 0, 1, are shown for 
Ma = as functions of q and cos#. Both drop quickly with increasing relative momentum 
between the two nucleons inside the deuteron. The wave function <^o(g, 9) shows a cosine-like 
behavior indicated by the straight line at q = connecting the highest point at 9 = with 
the lowest point at 9 = 180° through zero at 9 = 90°. This cosine-like behavior is confirmed 
to be true when considering the analytical angular behavior of the wave functions in Section 
1 1 [B[ The function <fii(q, 9), in contrast, displays a sine-like behavior, which also will prove 
to be its correct analytical form. It peaks at 9 — 90° and vanishes at 9 = and 180°. The 
maximum of (fn{q, 9) is larger than that of </??(g, 9). 

For Md = 1 and —1 the wave function components ip A d (q,9) are shown as functions of 
q and cos 9 in Fig. 0. They also drop quickly as the relative momentum between the two 
nucleons inside the deuteron increases. In the upper part of the figure we see that both 
(Po{q,9) and <^Q 1 (g,6 l ) vanish at 9 = and 180° but have a different sign for the other 9- 
values. At 9 = 90° fl(q, reaches its minimum whereas 1 (q, 9) reaches its maximum. In 
the lower part of the figure we see that <p\(q, 9) peaks at 9 = and vanishes at 9 = 180°. 
On the contrary ip^^q^) peaks at 9 = 180° and vanishes at 9 = 0. This angular behavior 
of (p\(q,9) and ip~^{q,9) suggests a relation between the two functions. This relation is 
explicitly given in Section pTTB| . For M d = 1 and -1 the maximum of ipf ld (q,9) is larger 
than that of (fo Id {q, 9). 

In Fig. ^ the deuteron densities p Md (q) as given in Eq. ( |2.13| ) for Md = (top row) 
and Md = 1 (bottom row) are shown. On the left side the two densities are displayed as 
functions of q and cos#, and on the right side as functions of the Cartesian projections of q, 
q x and q z . Since the wave functions are invariant under rotations around the z-axis, we show 
a vertical cut through the x-z-plane. The densities p°(q) and p 1 (q L ) have uniform angular 
distributions with respect to the azimuthal angle 9 for the momentum range shown, and 
thus the equidensity surfaces as function of the momentum between the two nucleons have 
a spherical shape. The densities are largest at small relative momentum q. Though not 
shown here the deuteron density for Md = —1 also has a similar shape. In all Figs. 1-3 the 
densities have their maximum at q = 0, indicating that the two nucleons being at rest with 
respect to each other is the most probable configuration for the deuteron. 

Using Eq. ( |2.26| ) we extract the usual S and D wave components. They agree very well 
with the ones obtained from standard partial wave calculations. 



III. FORMULATION II 



S 



A. Deuteron Wave Function in Operator Form 

In order to study the different spin orientations of the two nucleons in the deuteron in 
relation to the vector of relative momentum a representation of the deuteron wave function 
in operator form is ideal. It is also desirable to derive another set of two coupled, one- 
dimensional equations, in the basis of total helicity. 

In terms of the partial wave basis states given in Eq. ( 2.19Q the deuteron state has the 
well known form 

* d d ) = \t) E / d ^ 2 \l(tt)lM d ) Mq)- (3-1) 

l=0,2 J ° 

Here \t) indicates the isospin, which is for the deuteron. The explicit reference to t will 
be omitted in the following considerations. Again, we would like to point out that we work 
in the basis of total helicity, and thus our final expressions will differ from the ones given in 
Ref. ||. Carrying out the angular momentum expansion explicitly and using (q\ Ijj) = Yi^{q) 
one obtains 

V47T 



+ {|11) C(211; M d - 1, lM d )Y 2tMi -i(q) 
+ \10)C(211;M d 0M d )Y 2Md (q) 

+ |1 - 1) C(211; M d + 1, -1, M d )Y 2jMd+1 (q)} ^ 2 {q). (3.2) 
Inserting the explicit expressions for the Clebsch- Gordon coefficients [TT| leads to 

v^(q) = |lM d ) -4^Mq) 




2-M d )(3-M d ) 



20 



2-M d )(2 + M d ) v ... 

— *2M d yq) 
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(2 + M J(3 + MJ ywi( , nfeh) (33) 

Now we would like to express the wave function in a simple way such that 

\I/f«(q) = {c Mq) + c 2 Mq)} \lM d ) , (3.4) 

and where Co and c 2 are operators acting on the deuteron spin state | lM d ) . For this purpose 
we choose as example M d = 1 which yields 

*J(q) = |11) -!=Vo(<z) + I |H) \ ^Y x {q) - |10) \^Y 21 (q) + |1 - 1) \l~Y 22 (q) \ f 2 {q) 



4tt I V 10 y ' ' ' V 10 ' V 5 

1 MO)- (3-5) 
' zq" v ^ 
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ii) -j^Vt>(?) + {in) 0?o 2 + m-i) - [io> 3<mi + |i - 1) H) 



In the last step we expressed the spherical harmonic functions in terms of the spherical 
components of the momentum, g 1; q and g_i [|ll|]. Since the state \lM d ) in Eq. ( |3.4|) 
has already the correct transformation property under rotation of the deuteron state, the 
operators Cq and c 2 must be scalars under rotation. Those scalars have to be formed out of 
the spherical components of er(l) and er(2) which at the same time will connect the given 
states |1 — 1) , 1 10) and |11) to \lMd). Therefore we consider 

<r(l)-q<r(2).q|ll) = (a (l)q - ai(l)g_i - <r_i(l)gi) ' ' 



( < To(2)g -*T 1 (2)g_ 1 - -_i(2) gi ) 
g 2 |ll>-2g ogi |10> + 2g 2 |l-l>. 

1„2 



22 
1 1 

22 



The / = admixture can be projected out by subtracting -^q , which leads to 

'ff(l) • q<r(2) • q - ]-q 2 ) |11) = | (g 2 + q^) |11) - 2q ogi |10) + 2q\ |1 - 1) 



(3.6) 



(3.7) 



A comparison to the terms in Eq. ( |3.5| ) reveals that ^(q) can be written as 



47T 

Mq) + 



Mi) + 



<x(l)-q<7(2)-q--g 2 



<T(l)-q<7(2)-q--g 2 



Aq 2 V 2tt 
^2(?)} |H) 



Mq) I 11 ) 



where 



V'o(g) = 
$2(9) = 



1 



'47T 

3 1 



4g 2 V2T 



^o(g) 
^2(9) 



(3.8) 

(3.9) 
(3.10) 



We denote the expression in Eq. ( |3.8|) as 'operator form' of the deuteron wave function in 
momentum space. A corresponding expression in coordinate space can be found in Ref. J?|. 
In a fashion similar to the above derivation, one can show that the form given in Eq. (|3.8| ) is 
also valid for M d = and —1. Hence, the deuteron wave function in operator form is given 
in momentum space as 

1 



<T(l)-q<7(2)-q--g 2 



iP 2 (q)\\lM d ). 



(3.11) 



Here the positive parity is manifest, since \E^ /d (q) = ^ d (— q). It is a straightforward 
algebra to work out the normalization of ^f^f d ) as given in Eq. ( |3.11| ) and one obtains 



a 



dq(lM d \\Ml) + 

Mq) + 



<r(l)-q<x(2)-q--g 2 



<T(l)-q<7(2)-q--g 2 



^2(9) 

Mq)\\±M d ) 



An 
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(3.12) 



The last form is the standard normalization of the deuteron wave function in terms of partial 
wave components. In arriving at this result we used that 

<l«i| «r(l) ■ q CT (2) • q|!A« = { j _ ^ • «| I f (3.13) 

and 

/47T 
rfg (lM,| -(l)-q < r(2)-q|lM d ) = — q 2 . (3.14) 

As we shall show in Section [TV] based on the form given in Eq. (|3.11| ) one can express the 
angular dependencies of all possible spin orientations in the deuteron analytically. 



B. Analytical Angular Behavior of the Deuteron Wave Function and the Deuteron 

Eigenvalue Equation 

With the operator form Eq. (|3 . 1 1| ) at hand we revisit the deuteron wave function com- 
ponent in the helicity basis as given in Eq. 



^(q)= la (q;glA;0|^ 



(glA| «q| + <-q|) *f 
2<$lA|*f d (q) 

;i)- q< T(2)-q-^ 2 ]^ 2 (g)||lM d ). (3.15) 

The operator cr(l) ■ qer(2) ■ q can be expressed in terms of the total helicity S ■ q as 

<7(l)-q<T(2)-q=2(S-q) 2 -g 2 , (3.16) 
where S = | (c(l) + c(2)). Therefore, the helicity wave function component is given by 



2 {qlh\ |^o(g) + 



^(q) = 2(glAN^ (g) + 



2A 2 



2|Vo(g) + 

2 [Mq) + 
2 {Mo) + 



2(S • q) 2 - \q 2 
4 



$2(q)\\lM d ) 



q 2 Mq) (qlA\ lM d ) 



2A 2 -- 



2A 2 



(3.17) 



This shows that the angular behavior of the wave function component is given by 



JM d <hjl 



b M d A 



(6), where the d-matrix is explicitly given as [1 I 



d \l d A^) 



/ 1+cos 9 sin 9 1— cos 9 \ 

12 V2 2 \ 

sing n _ sing 

V2 s/2 

1— cos 9 sin 9 1+cos 9 

2 y/2 2 / 



(3.18) 
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Finally, we define an angle independent function $a(?) via 

4' 



$ A (g) = MQ) + 
1 



2A 2 

^o(g) + 



3 

3A 2 - 2 



9 ^2(?) 



1 



in 



=^2(9), 



(3.19) 



so that the deuteron wave function component can be expressed as 

#(q) = 2$ A (g)e^^^ A (0) 



(3.20) 



Here <£^ d (g, 5) are the wave function components which we determined previously numeri- 
cally and which are displayed in Figs. |l]-@. 

Employing the above given form of the deuteron wave function, we can derive a one- 
dimensional eigenvalue equation, starting from the eigenvalue equation for (p^ d (q,9) = 
y?A d (q)exp(-zM d 0). Inserting Eq. ( |3.20|) into Eq. (gT|) gives 

(£. - E d ) $ K{q)d\ Id M 

+ J dq' e -^-(^) {ivS°(q, c{)^ x {q')d l Mdl {9') + ^(q, q0$o(g')4, d o(^)} = 0, (3.21) 

an equation, which is valid for any direction 9. Choosing q = z simplifies the equation, since 
the azimuthal dependencies of the potential can be factored out as 



VfJ{qz, q') ee e iA ^V£p(q, q', 9'). 



(3.22) 



The d-matrix in the first term gives 5M d A, and the 4>' integration requires A to be equal Md, 
leading to 
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+7T jT dq' q' 2 £ d cos & {V™(q, q', ^(^i^) + \vS°(q, q', 9')^)d 



M d 0l 



= 0. 
(3.23) 



Choosing Md = 1 and leads to a closed system of two coupled equations for the amplitudes 
$1 and $o- The cos#' integration can be worked out independent of the amplitudes $a(<z')> 
so that Eq. fl3.23| ) is in fact a set of two coupled equations in one variable, namely q: 



Q_ 

m 



~E d ) $M d (?)+7T fd^^[vZMM^ + lv^MM^)} = ^ (3-24) 



with 



-1 



dcos9'V™ AI (q,q',9')d 



M d A'\ 



(3.25) 
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The set of two coupled eigenvalue equations Q3.24Q can be easily solved using the same 



method as described in Section IIP . The Gaussian grids for the q' integration and the cos#' 
integration in Eq. ( |3.25| ) are taken to be the same, and we obtain the same value for the 
deuteron binding energy, E d = 2.224 MeV. The solutions $o(<?) an d $i(<?) are displayed in 
Fig. |j. This figure shows that both functions are of the same magnitude for q = 0, and both 
drop by about one order of magnitude within q of ~ 200 MeV/c. $i(g) has its first node 
already for q « 300 MeV/c, while the first node of $0(9) occurs for q « 800 MeV/c. In 
general, the magnitude of $0(9) falls °ff slightly slower than the one for $i(g) as function 
of q. 

In Figs. ^ and ^| the wave function components (p^ d (q, 9) are obtained from numerically 
solving Eq. (|2~T8|) . With the help of Eqs. fl3~20D and (PH ) we can express their angular 



behavior as 

M d = : <p° (q, 9) = 2$ (g) cos 9 (3.26) 

<p\(q,0) = \/2$i(g) sin^ (3.27) 

M d = l:<p 1 {q,6) = -V2%(q)sm9 (3.28) 

^(g,0) = $i(g)(l + cos0). (3.29) 

Obviously, the angular behavior extracted numerically agrees with the analytical one. 

We mentioned in Section [II D| that the maximum of ip A d (q, 6) with A = M d is larger than 
that with A ^ M d . Eqs. ( gg§ )-( ggg ) show that the ratio \<PM d d (q, 9) max /ip^ Md (q,9) max \ is 
exactly \/2. This can be understood as follow. According to Eq. (|3.19|) the component 
ip^ d (q, 9) is determined for small q dominantly by ^o(q), i.e. the S wave. 

The analytical angular behavior of the deuteron densities given in Eq. ( p. 13 ) can now 
easily be derived. For M d = 0, 1 we find 

p°(q) = $l( q ) sin 2 9 + $ 2 (g) cos 2 6 (3.30) 
p X (q) = + cos 2 9) + Ul{q) sin 2 9. (3.31) 

From these expressions we can deduce that p°(q) and p x (q) are only perfect spheres for 
small q, where $o(?) and $i(g) are almost identical. For larger momenta the spheres are 
deformed according the ratio | $o( ( ?)/ ( ^ ) i (<?)!• 



C. Relation to the Conventional Partial Wave Representation 

Before completing the considerations about the analytic behavior of the angular be- 
havior of the deuteron wave function in the helicity basis, we want to make contact with 
the standard representation of the deuteron wave function. In Section |11 C| we derived the 



projection of the deuteron state on the partial wave basis. We ended up with Eq. ( |2.26|) 



and left the remaining quantum numbers j and I to be determined numerically. The wave 
function components (p\(q,9) together with their analytical angle behavior allows to calcu- 
late the projection and to determine the remaining conditions for their existence. Inserting 
Eq. (ET2QD into Eq. ( gg§ yields 
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Uq) = 2^(2/ + 1) 

f^dcose' \c{iij-Qiiy Mil ^ 

= ±5 jiy /ir(2l + l){c(llj; 011)$i(g) + ~C(Uj; 000)$ (g)} . (3.32) 
Here we use the orthogonality property of the d-matrix 

r\ 2 

j ^ d cos 9 d^ imi (9)d^ 2m2 (9) = - ^ S hj2 6^ 2 8 mim2 . (3.33) 

The projection exists only for a total angular momentum j = 1. Furthermore, the Clebsch- 
Gordon coefficients allow only I = and I = 2 and we obtain explicitly for the S and D 

wave 

M?) = jp/^{2$i(<z) + $o(g)} (3-34) 
ih(q) = -V^i^q) -$ (g)}, (3.35) 

which is consistent with Eq. ( |3.19| ). We extracted the S and D waves from Eqs. ( |3.34| ) and 
( |3.35| ) and found very good agreement with the ones obtained from a standard partial wave 
solution of the deuteron eigenvalue problem. 



IV. PROBABILITY DENSITIES FOR DIFFERENT SPIN CONFIGURATIONS 

The operator form of the deuteron wave function given in Eq. (|3.11| ) is an ideal tool 
to express probabilities for different spin configurations within the deuteron. This provides 
analytical insight into the shape of these configurations. As an example we choose a polarized 
deuteron with Md = 1. Cases of interest are if (1) both nucleons have their spins up, (2) 
both nucleons have their spins down, (3) one nucleon has spin up and the other has spin 
down, (4) one nucleon has spin up and the other has arbitrary spin orientation and (5) one 
nucleon has spin down and the other has arbitrary spin orientation. For these five cases the 
probability densities are given below. For clarity the final expressions are given in terms of 
the standard S and D waves. 

1. probability density for both nucleons having their spins up: 

p} T (q) = *J*(q)± [1 + (7,(1)] \ [1 + (7,(2)] *J(q) 

tPKi) + ^ (cos 2 tf - i) M<i)Mq) + 1 (cos 2 tf - £)% 2 2 (g)} • (4.1) 

2. probability density for both nucleons having their spins down: 
^(q) = ¥j*(q)i [1 - (7,(1)] \ [1 - a z {2)\ *J(q) 

° sin 4 ^ 2 2 (<?). (4.2) 



327T 
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3. probability density for one nucleon having spin up and the other having spin down: 
p T Vq) = *S*(q)i [1 + a z {l)\ ~ [1 - a z {2)\ ¥$(q) 

9 cos 2 6 sin 2 e^l(q). (4.3) 
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4. probability density for one nucleon having spin up and the other having arbitrary spin 
orientation: 

A4 (1) (q)=*i*(q)i[H-^(l)]«rJ(q) 
1 
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V>o(?) + 4| ( cos2 ^ - 3) Mq)Mq) 

j(cos 2 #-^ + cos 2 flsin 2 flj ^(g) 
= P ; T (q)+ P ; i (q). (4.4) 

5. probability density for one nucleon having spin down and the other having arbitrary 
spin orientation: 

^ a) (q)=*i*(q)i[l-a,(l)]*i(q) 
9 sin 2 ^ 2 (g) 



32tt 

= Pu(q) + Pu(q)- (4.5) 

In Figs. |5|-[5] those five different probability densities are shown. In each figure the left 
side displays the probability densities as functions of q and cos6*, whereas the right side 
depicts the probability densities as functions of q x and q z . The contour lines represent a 
vertical section in the x-z plane through equidensity surfaces. Rotating this section around 
the g 2 -axis gives a three-dimensional image of the equidensity surfaces. 

The probability densities for the first two cases, where both nucleons have the same spin 
orientations, are given in Fig. |5]. The top row represents pi*(q). The density peaks at q = 0, 
indicating that the largest densities occur at small momenta. This density has a spherical 
shape, since Eq. ( |4.1| ) is dominated by the S'-wave, and in the momentum range shown 
has little dependence on the angle 9. The figures in the bottom row represent pji (q). As 
Eq. ( (4.2|) suggests, this density is only determined by the deuteron D-wave times a function 
of the angle 6. Thus at q=0 it is zero, and reaches two maxima at \q ma x\ ~ 100 MeV/c 
along the g x -axis (9 = |). If a measurement could be carried out on a deuteron at rest the 
two nucleons would have momenta back to back perpendicular to the polarization axis of 
the deuteron. Rotating the vertical section given in Fig. |5|d around the z-axis will show a 
toroidal shape of the equidensity surfaces of the probability density in this configuration. 
For the image in Fig. 6, two equidensity surfaces, one with a high value, being closed in 
the section of Fig. 5d, and one with a small value are picked and rotated around the z-axis 
resulting in a torus, being cut open vertically. The surface of lower density is left half open 
at the outer side. The image displays a shape characteristic for the spherical harmonics with 
I = 2,m = 2. 
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For the case where the spins of the two nucleons point into opposite directions, the 
probability density is shown in Fig. [7[ According to Eq. ( |4.3|) , this density is also given 
solely by the deuteron D-w&ve and a function of the angle 9. It has four peaks of equal 
hight in each quadrant of the q x — g 2 -plane at \q x \ = \q z \ = q max cos(|). Rotating the vertical 
section in the x-z-plane around the z-axis will reveal a double toroidal structure. For the 
image in Fig. 8 two equidensity surfaces are picked and rotated around the z-axis, resulting 
in a double torus being cut open vertically. The inner tubes represent surfaces of higher 
density compared to the outer ones. The shape is characteristic for a spherical harmonics 
with I — 2,to= 1. Again, a measurement on the deuteron at rest would see in the maxima 
the two nucleons with momenta back to back pointing at 9 = 45°. 

For the remaining two cases given by Eqs. (j4.4|) and Q4.5|) , where only one of the two 
nucleons is polarized, the probability densities are presented in Fig. [|. The figures in the 
top row represent pj^^q). For the momentum range shown its properties are very similar 
to pj-f(q) given in the top row of Fig. |. The reason is that pjf (q) is larger than p^(q) and 
thus dominates. The figures in the bottom row depict pj^-^q). This density has the same 
maxima as pj^(q) given in Fig. [5]d, but a slightly different angular behavior. For a fixed 
q the changes with 9 are slower than for pj|(q). This is caused by the linear dependence 
on sin 2 9, whereas pjj^(q) has a quadratic one. A rendering of two equidensity surfaces is 
displayed in Fig. 10, which shows the different angular behavior when comparing to Fig. 6. 



V. SUMMARY 

As an object with internal structure it is tempting to investigate the deuteron properties 
three dimensionally. To that aim we study the deuteron properties in a representation based 
on the total helicity of the two-body system taken along the relative momentum of the two 
particles. Though originally developed for describing NN scattering, the method is general 
and can be used to solve bound state problems as well. 

We introduce deuteron wave function components in the helicity basis. They depend on 
the magnitude q of the relative momentum and the angle 9 of the relative momentum q to 
the z axis. Deriving an 'operator form' of the deuteron wave function one obtains insight 
into the analytical angular behavior of those components. 

We derived two sets of two coupled eigenvalue equations for deuteron wave functions. 
The first set of equations does not use any a priori knowledge of the quantum numbers of 
the deuteron, and the NN potential representation in helicity basis is used similarly as for 
NN scattering. As a consequence one has coupled two dimensional equations. In the second 
set of equations, the well known deuteron properties of a S and D waves content is build 
in and the analytical angular behavior of those amplitudes is taken into account explicitly. 
In this case, one arrives at one dimensional equations. This second derivation bears some 
similarity to partial wave methods. The two one-dimensional amplitude components are 
each linear combination of the standard 5* and D wave function components. 

The calculated deuteron binding energy determined in both ways agrees perfectly with 
the value determined in standard calculations based on partial wave representations of the 
deuteron eigenvalue equation. The newly defined helicity wave function components depend 
on Wigner's d-function and the q dependent part are linear combinations of the standard 
S and D waves of the deuteron. We display their properties for different projections of the 
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total angular momentum M^. Like for NN scattering we can connect the helicity amplitudes 
to the standard S and D waves and find perfect numerical agreement with the partial wave 
components determined in a standard manner. 

Finally we evaluate various spin and momentum dependent probabilities in a fashion, 
which is exact with respect to the angular dependence. This is made possible by using the 
'operator form' of the deuteron state. It is conceivable that in quasi elastic electrodesintegra- 
tion of the deuteron one may be able to see those momentum dependent spin distributions. 

Summarizing, we extended a recently introduced helicity representation for NN scattering 
to the NN bound state. This formulation leads to new forms for deuteron wave function 
components, which can be determined by two coupled equations. 
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FIGURES 

FIG. 1. The deuteron wave function components (fo(q,6) (a) and ipi(q,9) (b) in units 10 
MeV" 1 ' 5 as functions of q and cos 9. 



FIG. 2. The deuteron wave function components (p})(q,9) (a), ip\(q,9) (b), ip 1 (q,9) (c) and 
ifii 1 (q,9) (d) in units 10~ 3 MeV~ L5 as functions of q and cos 9. The relation between the different 



components is discussed in Section IIP 



FIG. 3. The deuteron density for = (a and b) and = 1 (c and d) in units 10~ 6 MeV -3 
as functions of q and cos 9 (a and c) and as functions of q x and q z (b and d). For the momentum 
range shown the deuteron densities display a uniform angular behavior. The contours represent 
equidensity lines along a vertical section in the x-z plane. 



FIG. 4. The wave functions |$o(?)| and |<E>i(g)| in units MeV 



FIG. 5. The probability densities pj-|-(q) in units 10~ 6 MeV -3 for both nucleons having their 
spins up (a and b) and pjj(q) in units 10~ 10 MeV~ 3 for both nucleons having their spins down (c 
and d). The contours represent equidensity lines along a vertical section in the x-z plane. 



FIG. 6. Two selected equidensity surfaces of the probability density pjj(q) for both nucleons 
having their spins down. The image is created by rotating two of the equidensity lines of Fig. 5 d 
around the z-axis. Note that the z-axis is stretched with respect to the other two axes. 



FIG. 7. The probability density pjj(q) in units 10 -10 MeV -3 for one nucleon having spin up 
and the other having spin down. The contours represent equidensity lines along a vertical section 
in the x-z plane. 

FIG. 8. Two selected equidensity surfaces of the probability density pjj(q) for one nucleon 
having spin up and the other having spin down. The image is created by rotating two of the 
equidensity lines of Fig. 7 b around the z-axis. 

FIG. 9. The probability densities pj^(q) in units 10~ 6 MeV~ 3 for one nucleon having spin up 
whereas the other having arbitrary spin orientation (a and b) and pj^^q) in units 10 -10 MeV -3 
for one nucleon having spin down whereas the other having arbitrary spin orientation (c and d). 
The contours represent equidensity lines along a vertical section in the x-z plane. 
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